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Scattering of color dipoles: from low to high energies 
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A dipole-dipole scattering amplitude is calculated exactly in the first two orders of perturbation 
theory. This amplitude is an analytic function of the relative energy and the dipoles' sizes. The 
cross section of the dipole-dipole scattering approaches the high-energy BFKL asymptotics starting 
from a relatively large rapidity ~ 5. 
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I. INTRODUCTION 

It is well known that the effective degrees of freedom 
at high energies are so-called Wilson lines - infinite gauge 
factors ordered along the velocities of the colliding par- 
ticles (for a review, see Ref. |l|). Let us consider the 
classical example of the scattering of virtual photons in 
QCD. The photons decompose into quark-antiquark pairs 
which interact by exchanging gluons. At high energies, 
quarks move very fast so that their propagation in the 
background fields of the exchange gluons is reduced to 
the gauge factor of U^{x±) 



— pig I^^dun^A''(un+x±) 



(1) 



ordered along the classical trajectory of the particle, i.e. 
a straight line coUinear to the velocity n. Here xj_ is 
the transverse position (an impact parameter) of the fast 
quark, which does not change in the collision. The propa- 
gation of a quark-antiquark pair is described by the color 
dipole (3] formed from the two Wilson lines 



W^"(a:i,2/J-) = Tr{C/"(xj.)f/"t(y^)} 



(2) 



where x± and y± are the transverse positions (impact 
parameters) of quark and antiquark. The high-energy 
7*7* scattering reduces then to the dipole-dipole ampli- 
tude integrated over dipoles' sizes and separations (see 
Fig. P): 

A{s) ^ d?a^(fb^ I^{a^)l''{h^)T{a^,br,s) (3) 

where I^{a±) is the so-called "impact factor" (see Ap- 
pendix A) and T(a±, b±; s) is the dipole-dipole scattering 
amplitude: 



T{a^,b^;s) = 

-% f d2zi(W^"(aj. + zj_,zj_) W%b±,0±)) 



(4) 
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FIG. 1: High-energy scattering as a dipole-dipole amplitude. 
Wavy lines are gluons and dashed lines are Wilson gauge fac- 
tors Q 



Here e is the unit vector in the direction of the motion of 
the second virtual photon. The collision energy is related 
to the relative rapidity rj (= angle between the dipoles) 
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= n ■ e = cosh 77 



(5) 
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The energy dependence of the 7*7* amplitude is governed 
by the evolution of the dipoles with respect to the slope 
of Wilson lines determined by 77 "7, 'H^ . 

We see that the dipole-dipole amplitude is an "elemen- 
tary" high-energy scattering process. Moreover, in theo- 
ries without asymptotic particle states (like N=4 SYM) 
the dipole-dipole scattering is the only way to access the 
high-energy behavior of amplitudes. Also, there have 
been many attempts to estimate the high-energy behav- 
ior of hadron-hadron amplitudes using non-perturbative 
models for the dipole-dipole scattering (see Ref. for a 
review). For hadrons, the impact factors cannot be cal- 
culated in pQCD; however, they can be related to the 
hadron wave functions. 

The usual phcnomenological approach to dipole-dipole 
scattering in pQCD is to take the light-like dipoles. How- 
ever, matrix elements of the light-like Wilson operators 
contain longitudinal divergencies. In the LLA we cut 
these divergencies "by hand" |3, . This prescription - 
take the light-like 
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dipoles and cut the divergencies by hand - seems to 
avoid ambiguities in the next-to- leading order, tooQ- 
However, to go beyond perturbation theory one has to 
consider the scattering of the ofF-hght-cone dipoles. The 
dipole-dipole amplitude is then a function of the dipole 
separations and the angle 9 (~ relative rapidity "q) be- 
tween the dipoles. If this amplitude is an analytic func- 
tion of 6, one can calculate the dipole-dipole scattering 
(a correlation function of two infinite rectangular Wil- 
son loops) in the Euclidean region and continue to the 
Minkowski space. In the recent series of papers [Tol [Tll | 
the dipole-dipole scattering was calculated in the Eu- 
clidean space using the AdS-CFT correspondence and 
continued analytically as a function of the angle 9 to 
the Minkowski space, where 9 — *■ irj. Another example 
of such analytic continuation is the calculation of the 
instanton-induced dipole-dipole ampltitude [T^ . 

To be able to obtain the high-energy behavior from 
the Euclidean amplitudes it is crucial to verify that the 
dipole-dipole amplitude is an analytic function of 9. In 
Ref. this statement was proved for the correlation 
function of two Wilson rectangles with a finite longitudi- 
nal length T (see also Ref. Il^ . At finiteT, the amplitude 
is an analytic function of the angle 9 between the rectan- 
gles. In the limit T oo, one recovers the dipole-dipole 
scattering amplitude. However, whether the T — > cxd 
limit commutes with the analytic continutation from the 
Euclidean to Minkowski space is an open question. We 
demonstrate with an explicit calculation in the first two 
orders of perturbation theory that the dipole-dipole am- 
plitude is, indeed an analytic function of the angle be- 
tween the dipoles (as well as of the dipoles' sizes). This 
is the first important result of our paper. We see that 
the analyticity survives the T ^ oo limit in the first two 
orders of perturbation theory, and presumably it hap- 
pens in all higher orders as well so that the dipole-dipole 
amplitude appears to be an analytic function of 9. 

The second result of our paper concerns the rate at 
which the dipole-dipole amplitude approaches the BFKL 
asymptotics. At sufficiently high energies, the cross 
section for the "unpolarized" dipole-dipole scattering is 
given by the BFKL formula 

a(a,6;77)~^ J rf^(^) (6) 

where x(j^) = — 9fti^(l/2 + iv) — C is the position of the 
"hard pomeron" (for a review, see Ref. [T6t) . The 
asymptotics of the 7*7* scattering is then given by Eq. 

The formula © is obtained in the so-called lead- 
ing log approximation (LLA), when as ^ l,Q;slns 
asTj ^ 1. The main problem with the LLA result ® 
is that its power behavior s'^'^'^"^'^'^ violates the Frois- 
sart bound In^ s at asymptotically large s. Thus, the 
BFKL pomeron gives only the pre-asymptotic behavior 
at "moderately high" energies. It is a common belief 
that the true asymptotics at s ^ 00 (where alns 3> 1) 
comes from the unitarization of the BFKL pomeron. De- 
spite numerous efforts, the high-energy asymptotics sat- 



isfying the Froissart bound is still an unsolved problem 
in QCD (for a recent discussion see Refs. 0, 0, IT^ . 
In this paper, we address a different problem (important 
for mid-energy accelerators) - when does the LLA asymp- 
totics start making sense - at 1, 10 , or 100 GeV? To get a 
complete answer to this question one should calculate the 
amplitude exactly and compare it to the BFKL asymp- 
totics Since it seems to be an impossible task at 
present, we calculate the dipole-dipole scattering in the 
first two orders of perturbation theory exactly and com- 
pare it to the asymptotic form. Our main result is that 
the asymptotics starts rather late, at 77 ~ 5, which trans- 
lates into y/s ~=10 GeV for the scattering of dipoles with 
the p-meson size a ^ O.S/m. 

The paper is organized as follows. In Sect. 2 we calcu- 
late the dipole-dipole amplitude and the cross section in 
the leading order of perturbation theory. In the second 
order, we have calculated independently both the dipole- 
dipole amplitude and its imaginary part proportional to 
the cross section. Since the structure of the diagrams 
and the result is much more transparent for the cross 
section, in Sect. 3 we present only the results for the 
amplitude in the second order and leave the details un- 
til Sect. 5, where we give the detailed calculation of the 
dipole-dipole cross section. We discuss the effects due to 
the running coupling constant in Sect 4. Sect 6 contains 
the numerical estimates of the dipole-dipole cross section 
and Sect. 7 is devoted to the discussion of our results 
and an outlook. The explicit form of the second-order 
amplitudes is given in the Appendix. 

II. DIPOLE-DIPOLE SCATTERING IN THE 
BORN APPROXIMATION 

We consider the dipole-dipole scattering amplitude de- 
fined as follows 

T{xi±,X2±;x3j_,X4,±;r]) (7) 

= -l(VK"(xii,.X2i)I^^(x3±,4i)) . 

Strictly speaking, the Wilson lines forming the dipole are 
connected with gauge links at infinity, so 

T{xi±, X2±; X3±, X4±;t]) 

= -tlim {W^{xi^,X2±)W^{x3±,Xii_)) (8) 

T — >oo 

where is a rectangular Wilson loop of longitudinal 
length T 

T T 

W:^{xi^,X2i_) = Tr{U^{xii_)[xii_ - —n,X2i_ - —n] 

X U'Ji-\x2i_)[xii_ + -n,X2i_ + -n]} (9) 

and [/^(xij^) = [^n + xi^, — '^n + xnJ\ etc. Hereafter, 
we use the notation 

[x,v\ = Pe^f^o '^"(^-!/)''^M("^+(i-")2/)sec : (10) 
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for the straight-line ordered gauge hnk suspended be- 
tween the points x and y. 

Since the gluons reduce to the pure gauge fields at in- 
finity, the precise form of a contour connecting the points 
xi±_ ± and X21. ± does not matter. Moreover, we 
use Feynman gauge for the gluon propagator, and in this 
gauge the links at infinity do not contribute to the am- 
plitude. Thus, we omit them in what follows to simplify 
the notations. 



A. Cross section of the dipole-dipole scattering 
and optical theorem for dipoles 

The total cross section of the scattering of the dipole 
of size a on the dipole of size h may be defined as 

CT(a_L,5_L;77) = J dzj_<,{a± + z±] z±,b±]r-i) (11) 

where 

,{x,±,v) = ^(0|W^';^(a:ix,a;2±)W^,^/(a;3±,a;4i)|X) 
X {X\W;i{xi^,X2±)Wf,,{x3^,x^^)\0) (12) 

is the total amplitude of the dipole-dipole transition into 
hadrons. Here the summation goes over all the interme- 
diate states except for the vacuum and 

W^j{xi^,X2i_) (13) 
= lim iU^{xi±)[xi± - —n,X2± ~ —n]UJj^\x2±)) ■ 

T^oo \ Z Z / ij 

If we include the intermediate vacuum state in the r.h.s. 
of Eq. H12I) , we get due to the completeness relation. 
Separating the non-interacting contribution in a usual 

way 

{0\W{xi^,X2±)W''{x3±,Xi^)\0) - + iT{x,, ri) 

(14) 

we get the "optical theorem" for the dipole-dipole scat- 
tering: 

3T(xi_L,a;2i;a;3_L,X4_L;ry) = ^<;(xij^, a:2±; 2:3^, a;4±; r?) . 

(15) 

It is worth noting that the high-energy 7*7* cross sec- 
tion reduces to the dipole cross section (fl^ integrated 
with impact factors, similar to the Eq. (j^J for the ampli- 
tude: 

Ms) = \ I d^a^d^b^ I^{a^)I^{b^)a{a^,br,v) (16) 

We will calculate the dipole-dipole cross section in two 
ways: directly as the r.h.s. of Eq. H12|) or via the optical 
theorem as the imaginary part of eq. Q. 
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FIG. 2: Dipole-dipole scattering in the lowest order of per- 
turbation theory 



B. Lowest order dipole-dipole amplitude 

In the leading order of perturbation theory, the correla- 
tion function of two Wilson-line operators is proportional 
to the (massless) two-dimensional propagator: 

{U-ix^Wiy^) = -{U-ix^)U^Hy±)) 

= ig^ cothrj f^lh^ e'C^'^-y^^ (17) 

where (a, b)± denotes the (positive) scalar product of 
two-dimensional vectors a± and b±. (For brevity, we 
use the ?i-inspired notations (f"fc = ^^^j^^ and (f ("^(fc) = 

(27r)"(5(")(fc)). 

The relevant Feynman diagrams are shown in Fig. |21 
The contribution of these diargams has the form: 
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-i(ki ,X2) ± \( i(ki ,X3) ± 



(18) 



X ^gi(fc2,a;i)i _ gi(fe2,a;2)i'j^g-i(fe2,2:3)i _ g-i(fe2,2;4)i "J 

The integral over the transverse momenta can be per- 
formed exphcitly, resulting in 



T{xij_, X2±; X3^, x^±^; r]) 



2 /j^-L^L' ' 



9 9 
X23X]^4 



where X12 = xi± — X2± etc. (cf. Ref. il9^) 



C. Dipole-dipole cross section in the Born 
approximation 

As we mentioned above, the dipole-dipole cross section 
can be obtained in two ways: via the optical theorem as 
the imaginary part of Eq. p9l) or, directly, as the r.h.s. 
of Eq. (|14|l . Since the amplitude in the lowest order in 
as is purely imaginary, the "unintegrated cross section" 
(ini is 



1 



coth^, (ln±l^) (20) 

•''23-^14 



For future use, it is instructive to calculate 'i{xi±,r]) di- 
rectly as the r.h.s. of Eq. p2|) . In this case, Feynman 
rules for the calculation of the cross sections can be re- 
produced by a functional integral over the double set of 
fields (see e.g. |20|): (+) to the right of the cut and (-) 
to the left with the propagators 



Al^{x)A'^{y) 
A'l^(a;X(y) 



dk 
dk 



■Jj.u:ab 
ik(x-y) " 



167r4 



(21) 

e-'H^-y)gt'''S''^2TTS{k^)e{ko) 



These correspond to the Cutkovsky rules for the cross 
sections. 

In these notations, the dipole-dipole cross section reads 

(cf. 113) 

<;{x,r,v) = (Tr{W^("_)(xii,X2i)W^(';)(xii,X2±)} 

X Tr{VK('L)(a;3±,X4i)W^f+)(a;3±,X4i)}> (22) 

The relevant diagram is shown in Fig. Q. To emphasize 
that here we have only the gauge links at T = —00 we 
draw them explicitly. 

The leading-order correlation function of Wilson lines 
in (-I-) sector is given by Eq. p7|) and in {—) sector by 
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FIG. 3: Cross section of the dipole-dipole scattering in the 
lowest order of perturbation theory 



Eq. (|17|1 with a different sign 

(x^)t/f+) (2/^) - -{Ul\){x^)Ull^{y^)) 



ig coth?7 



d^k 



in^kl 

Using Eq. H23|l it is easy to see that 
4 



(23) 



^2 _ ^^{xi±i X2±--, X3±, Xij_; v) 



(f%jjfk2^, 

2 l2 \^ 



2-i^(^-i{ki,xi)± _ g-i(fei,a;2)j 



i{ki,X3}± _ gi(fci ,2:4)^ -J ^gi('£2,a;i)j 

iik2,X3)± _ —i{k-2,X4,)± 



^i{k2,X2)j 



(24) 

where the 2* = 16 terms correspond to 16 diagrams in 
Fig. 131). Integrating over k±, we reproduce the optical- 
theorem result (I20|l . 



III. SECOND-ORDER AMPLITUDE 

In the next order of perturbation theory there are too 
many diagrams to present them all - rougly speaking, we 
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FIG. 4: Typical 
amplitude. 



'connected " diagrams for the dipole-dipole 



must take the diagrams in Fig. |21and add an extra gluon 
line in all possible ways. Let us consider the diagram in 
Fig. 12^. If we insert the gluon line in all possible ways 
to connect the left and the right parts of the diagram in 
Fig. we get the diagrams shown in Fig. ^ 

Also, we can insert a gluon line in the left or right part 
only and get the "disconnected" diagrams shown in Fig. 

ini 

Not all the diagrams for the correlator of two Wilson 
loops O contribute to the dipole-dipole scattering am- 
plitude. There are two exceptional classes of diagrams 
proportional to the total time of the evolution T. First, 
there are those shown in Fig. Elq-t, which describe "mass 
renormalization" of Wilson lines. Their contribution has 
the usual form of SmL, where Sm is the self-energy cor- 
rection and L ~ 2T is the perimeter of the Wilson loop 
corresponding to the dipole. Similarly, the diagrams in 
Fig. 21 m-p contribute to the "binding energy" e of the 
dipole and lead to eT terms. These two contributions 



exponentiate so that one obtains: 

-i{TT{U^{xi)U:^Hx2)} TT{U^ixs)U^Hxi)}) (25) 



-4iSmT -i{ei2+<:34,)T 



T{xi±, X2±,X3±,X4:±, m ^) 

where T(xn_, rj; T) describes the dipole-dipole scattering 
at the time T. The resulting amplitude 



T{xi±,ri)= lim T{xi±,T];T) 



(26) 



= lim -i{TT{Unxi)U:^Hx2)} TvIUUx^WtHxi)} 

^ g4i(5mT+i(ci2+e34)T 

is finite as T ^ oo. In the second order of perturbation 
theory, the multiplication in r.h.s. of eq. (|26() reduces to 
the subtraction of the corresponding terms (ASmT and 
2eT). The result is that the diagrams in Fig. 2|m-p and 
Fig. |S1 q-t are left out and in the diagrams in Fig. ^ 
i-1 and Fig. |5lm-p the color factor —jpj- is replaced by 

Thus, we are left with the diagrams shown in Fig. 
j and [SJ^-r. Also, there are ~ 20 times more diagrams 
which are obtained from the addition of an extra gluon 
to the graphs in Fig. 2b,c,...s. The calculation of these 
diagrams is standard but lengthy. In Sect. we will 
present some details of the (independent) calculation of 
the imaginary part of this amplitude where the structure 
of diagrams is much more transparent. Here we will only 
give the final result for the amplitude and discuss which 
classes of diagrams contribute to it. 

We have performed calculation of the dipole-dipole am- 
plitude in both the Euclidean and Minkowski spaces. 
In the Euclidean space, this amplitude is a correlation 
function of two Wilson rectangles at angle 6 such as 
cos9 ^ n ■ e. In both cases, the amplitude has the 
form (in Euclidean space rj = i9): 



T{X^±,1^) = 



1) 



An 



coth rj 



(27) 



X (e" 



-i{k[,X3) 



i{k[,X4)-^^^i(k2,xi) _ ^i{k2,X2)-^ 



X (e^(fc2,^3)_L _ ^^{k'2,X,)^^ ) ^2(^^ ^ ^/ 



A(k,;r]) (N^ - i)cothr] 

' ITT 



(fci + fci)6 ' N^(ki + k[)l 

+ t^{ki+k[)&^{k2 + k'2) 
k\ 



"-l"- 2 



1^2'^ 1 



in in — 4 

y ^ 



-J3)2 p'^{k2-p)' 



X coth 77 / d'^p ( ^ 
j Vp^iki 



Here /i is the normalization point of the MS scheme. The 
explicit form of the function j4(fci, k[,k2,k'2] rj) is given in 
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- "1 - >^2 



^ X-j ^ 



(a) 



(d) 



(b) 



(c) 



(f) 



Section). It is easy to see that the apparent divergency 
at fc + fc' ^ in the gluon reggeizatoon terms cancels 
with the corresponding contribution to A{ki; 77) (see Eq. 
EUl- 

The main conclusion from the Eq. (|27|l is that the 
dipole-dipole amplitude is an analytic function of the an- 
gle T]. It is easy to check that the functions A become 
real as r] ^ i9 (see the Appendix B) and so does the 
amplitude iT(xi±, iO). 

It is worth noting that in the LLA (as rj 00) we 
reproduce the first iteration of the BFKL kernel 



(g) 
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Cm) 



(p) 



(s) 
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FIG. 5; "Disconnected" diagrams for the dipole-dipole am- 
plitude. 



the Appendix B. Note that they are analytical functions 
of 77 and kfj_ with singularities at 77 — > corresponding 
to the bound state of two parallel dipoles. 

The r.h.s. of Eq. H27I) consists of two terms. The first 
term comes from the diagrams of the type shown in Fig. 
^ (The — 4 term is the odderon contribution com- 
ing from the structure d°''"^d°''"^). The first expression of 
the second term comes from the "gluon reggeization" di- 
agrams in Fig. Ela-d, while the second one (~ f — |^) 
comes from the gluon and quark self-energy shown in Fig . 
|S1 e-h. The diagrams in Fig|Sl i-1 vanish (cf. Ref. 123 1. 
while the terms coming from Fig[Slm-p are the pure di- 
vergency ~ J which also vanishes in the framework 
of dimensional regularization(see the discussion in next 



T(xij_,77) 



^^00 ig^v N,{Nl - 1) 
47r 4 



(28) 



(fcl - k^f 



t^{ki+k[-k2-k'2) 



ki k2 



^2 ^1 



1 2 1 2 



fci"(fci +fci)2 



IV. RUNNING COUPLING CONSTANT 

As we shall see below, the diagrams contributing to 
the renormalization of coupling constant are somewhat 
unusual: as in the case of a heavy-quark potential , 
the coefficient — in front of In/i comes from 

both the UV- and the IR-divergent diagrams (cf. ^22]).To 
see how it happens, it is instructive to start with the 
scattering of the light-like dipoles where the diagrams 
for the renormalization of coupling constant are the usual 
UV-divergent vertices and self-energies shown in Fig. 

As we mentioned above, the scattering amplitude for 
the light-like dipoles is not a well-defined quantity be- 
cause of the longitudinal divergencies corresponding to 
rj = 00 coming from the diagrams in Fig. 2|a-h and Fig. 
O a-d. If we cut this divergency "by hand" , we get the 
asymptotical contribution, i.e. Eq. (|28|l . Apart from 
this contribution, the only non-zero diagrams are those 
shown in Fig. [Sje-h. They lead to the renormalization of 
coupling constant in Eq. p9|l . 

Standard computation of these diagrams yields 



V3 3 



2ni 



2Nr 



9 



167r 



,ln^ (29) 



where ( |7V< 



^|7Vc — —^j comes from the diagram in Fig. 

and the similar diagram with the quark loop (see Fig. 
Od), while 2A'c comes from the diagram in Fig. (Sja. It 
turns out that the diagram in Fig. E); vanishes, while the 
diagram in Fig. Eli is irrelevant since it contributes only 
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(b) 



(c) (d) 

FIG. 6: Running coupling constant for the light-like dipoles 




(a) 



(b) 



FIG. 7: IR divergent diagrams contributing to the renormal- 
ization of the coupling constant 



to the "mass renormalization" 5m of the Wilson line (see 

Eq. (uni). 

Surprisingly, the structure of the diagrams contribut- 
ing to the running coupling constant is quite different for 
even slightly off-the-light-cone dipoles. In this case the 
diagram in Fig. is the same as for the light-like dipoles, 
the diagram in Fig. vanishes, and the diagram in Fig. 
EJ; is a pure divergency of the type / , which must be 
set to zero in the dimensional regularization approach. 
Following Ref. it is convenient to write down this 
term as 



87r2 



i-d d-4: 



(30) 



The UV pole together with the pole coming from 

the diagram in Fig. Eji forms 4T2m^(^'T^'^ ~ l*^/)' 
while the IR pole cancels with similar poles in the 
diagrams in Fig. [7| 

Due to "cancellation" between the UV and IR contri- 
butions to Fig. in order to reproduce the running- 

Nc — |rt/j in front of Infj,^ in 



couphng coefficient 



the dimensional regularization, one should take into ac- 
count not only the "usual" UV divergent diagrams in 
Fig. I^lbut also the IR-divergent diagrams in Fig. [T] (cf. 

[2^). As a result the coefficient ^|7Vc — in front 

of In /i^ comes from the diagram in Fig. while 



^gj^" In n"^ comes from the diagrams in Fig. [7| (and all 
other diagrams in Fig. vanish). 

This mixture of the IR and UV divergencies could have 
been a source of a potential problem. The cornerstone of 
the BFKL approach is the assumption that all the Ins 
come from the longitudinal integrations, while the inte- 
grations over the transverse momenta are convergent at 



where m is of order of masses (or virtuali- 



ties) of the colliding particles. (After the summation of 
the BFKL ladder we have "diffusion" in the transverse 
momenta leading to In fc^ ~ ^In but such terms do 
not spoil the "power counting" in (Ins)"). This assump- 
tion is not true for the diagrams leading to the renor- 
malization of the coupling constant since they are UV 
divergent. However, the structure of these diagrams is 
very simple - they are 1-loop vertex and self-energy cor- 
rections, and after the subtraction of the poles the 
integrals over fc^ are bound with /i^. In our case, the 
structure of the integrals over the transverse momenta is 
more complicated due to the mixture of the IR and UV 
divergencies. There are individual diagrams where the 
upper cutoff in the integrals over the k± is the energy 
^/s rather than - or the normalization scale 11. Such 
contributions would give the additional non-BFKL Ins 
terms coming from the logarithmical integrals over the 
transverse momenta J . Fortunately, we will see be- 
low that such terms cancel in the sum of all diagrams and 
the remaining transverse integrals are cut by the dipoles' 
sizes (or by the normalization point fi). 



V. CROSS SECTION 

A. Calculation of Feynman diagrams 

In the case of the cross section H12(l . the structure of 
the diagrams is more transparent and we will present 
some details of this calculation. In this order, there are 
two types of Feynman diagrams: one with gluon emission 
and one without. Typical diagrams of the first type are 
shown in Fig. |H1 The sum of these diagrams is propor- 
tional to the square of the finite-energy Lipatov vertex 
i^(fci, k[)Lfj_{k2, k2). The Lipatov vertex 



L{ki,k[) = (kl ~ k[){ne) + \2{k[n) + -^^k'l 



(fcie) 



(k[n) 



(31) 



is a sum of the emissions of the real gluon with momen- 
tum (kl + k'l) from the left part of the diagrams in Fig. 
13 in all possible ways. This vertex is gauge-invariant: 



{ki+k[)^'L{ki,k[)^0 



(32) 



At infinite energies (ne) —^ 00 the expression (|31|l reduces 
to the usual asymptotic form of Lipatov vertex (see e.g. 
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FIG. 8: Typical diagrams for the dipole-dipole cross section 
with the emission of a real gluon 



The sum of the diagrams in Fig. |H1 has the form 
4 



^ d'^kid'^k[d'^k2d*k'2 



(33) 

^S'^^^ki + k[ - k2 - k'^) 



USn'^kfk'lklk 
X d{ki + k[fe{ki + fc;)o 6{ki ■ n)S{k2 ■ n) 

X |'gi(fcl,a:3)i 



j,i(fcl,X4)j 

e-'('=='"^)^)[-L(fci,A:l).i(fc2,A:^) • 
The exphcit form of the product of two Lipatov vertices 



'L{ki,k[) ■ L(k2,k'2) 



(34) 



4(fc'L + fcL) - 2(ne) [^(fc'? + fc'^) + ^(fc^ 



{nef 



12 12 



'^le 



where fcie = (fci • e) and k[j^ = (fcj^ • n). Apart from the 
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FIG. 9: Typical diagrams with an extra virtual gluon 



diagrams with the emission of a real photon, there are di- 
agrams with an extra virtual photon. If, for example, we 
take the diagram in Fig. and insert an extra gluon line 
in the left (— ) sector in all possible ways, we get the dia- 
grams shown in Fig. We do not display the diagrams 
corresponding to 6vn or (5e, see the discussion in Sect. El 
The contribution of these diagrams vanishes here anyway 
since the phase factors e~*(2™+^)^ and e*(^™+'^)^ coming 
from the left and from the right of the cut cancel out in 
the cross section. Again, this leads to the replacement 

of the color factors ~2F c_f = — in the 

diagrams of the Fig|SJzi and Fig^-i type. 
The calculation yields: 

^2_ j^ '?(3^iJ-:^2±;a;3±,a::4^;??)FigEl 
/A^cCothS/ ( d^kx^d^k2^ 



(35) 



47r 



I67r4fc2^fc2^ 

i(fei,X2)i-jj-gi(fei,2:3)j 



X (^g«(*:2,2;i)i _ gi(fe2,2:2)i ■jj^g-i('£2,2;3)j 



^i(fei,a;4) 

-i(k2,X4)j 



) 



1] coth f] 



fr2 



k"^ 



-Ifl 

47r V 



1 

4^V kj^ 



4^2p^ V(fci_p)±2 ' (fci-p)±^ 
2 



In 



'^2J 



2nf_\ 
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where the first "gluon reggeization" term comes from the 
diagrams in Fig. |21a,b and the second from Fig. |51n,o. 
As in the case of the amplitude, the diagrams in Fig. 
j-m vanish and the diagrams in Fig. |^ h,i give rise to a 
pure divergency ^ / that does not contribute in the 

framework of the dimensional regularization. (Instead, 
it leads to the "mixing" of the IR and UV singularities, 
see the discussion in Sect. II Vp . Finally, the diagrams in 
Fig. l^-e cancel with the corresponding diagrams with 
an extra gluon to the right of the cut. 

It is worth noting that the diagrams in Fig. inif,g should 
be regularized by finite T before calculation. After sum- 
mation of all the relevant diagrams (e.g. Fig. |^f,h and 
Fig. |S|g,m,n), the limit T ^ oo becomes regular. 

Performing the remaining integrations over longitudi- 
nal momenta in Eq. and the integration over p in Eq. 
1351 one obtains 



167rtanh 77 



(t'''kict^k2d'^k[ct^k'2 (36) 



X (^^i(k2,xi) _ ^i{k2,X2)-j^^i{k'2,X3)± _ ^i(k'2,X4)±-j 



X S-^{ki + k[-k2-k'^] 



" iki + k[)^ 



167rtanh j] 



X (e~*('"'i'^i) — ^-iiki.X2)^^^i{ki,X3) _ gi(fci,a;4)-j 

i{k2,Xi)\ 



X (^g*(fe2,2:i) _ ^i{k2,X2)-^^^^-i{k2,xz) _ g-i(fc2 ,0:4) -j 
kl 



^2 



p)^ p^{k2—p)' 



1 



/ 5 Uf 



1 



Here 



(37) 



where oi, a2, 03, 04 and 05 correspond to the first, second, 
third, fourth and fifth terms in Eq. I|34() . respectively. It 
should be mentioned that the fifth term in the r.h.s. of 
Eq. (|34ll is IR divergent. After regularization, in addition 
to ar-,[ki; rj) it produces the In fi^xi2X34 term which forms 
— together with the contribution from the self- 
energy diagrams in Fig. |3n,o. 

The explicit form of the functions involved is: 

ai{ki,k2,k[,k'2;r]) = ^Ai{ki,k2,k[,k!2;r]), (38) 

where Ai can be found in the Appendix B. We see that 
the "unintegrated cross section" lp!Hjl is given by the 
imaginary part of the amplitude (|27|l (see the optical the- 
orem for dipoles (|15|) 'l. 



The cross section of dipole-dipole scattering is given 
by the integral 



cr(a, 6; 77) 



(fiz (^{a + z, z, b, 0; rj) 



(39) 



B. Asymptotics of the cross section 



As r/ — > 00, the cross section (|39|) reduces to 



<Tasy(a, &;s) = 



(40) 



'd^kcf^k' 

fc2 fc'2 



. • 2 (fc>a) .2 {k',b) 
4 sm sm 



1 



2tt 



fc2 



p^{k — py 



1 , /5 2nf 
— In ^ 

47r fc2V3 3A^c 



ln^(^-^) + ^\nfi'^ab[s-{k + k') 



21] 



\n{k + k'f/e (2e 



2-K 1 (k + k'f (k + k')^~k^ 



\n{k + k'Y/k '^ f2k' 
{k + k')^-k 



e + k'^ 
k^k'^ 



+ 0(6-") 



It is easy to see that the integrals over the transverse 
momenta converge at fc^ ~ a~^. In the individual dia- 
grams, the integrals over fc^^ diverge (which means that 
they would be cut by s in the corresponding exact expres- 
sions (|39|l and (|36() ) , but in the sum of all diagrams such 
terms cancel (see the discussion in Sect. II V|) . The struc- 
ture of the asymptotic result Eq. is 'qa{a, b)+a'{a, b). 
The term ~ 77 is the first iteration of the BFKL kernel. 
After integration over k and k' , it gives the first term 
in the expansion of r/"='<:('^) in powers of as- The second 
term may be called the "dipole impact factor" . Indeed, 
the standard representation of the BFKL asymptotics of 
the cross section has the form 



'asy 



I{a,k±)T]^^^^^I{b,k'^) (41) 



where /^bfkl is the BFKL kernel (the LLA kernel + 
NLO BFKL kernel +...). The kernel terms come from 
the central region of the rapidity, while the impact factors 
come from the fringes of the longitudinal integral where 
the momentum of emitted gluon is almost coUinear to n 
or e. In this paper, we do not access the NLO BFKL 
terms (those would correspond to contributions to 



the amplitude 



In s) so that all the non-LLA terms 



in Eq. (|4U|1 should be included in "impact factors" of the 
dipoles, namely. 



I{b, fci) = sin^ + 2ask^ 



d^k' 



ln(fcH-fc')Vfc2 /2fc2 
{k + k')^ - k^ [~k 



2 



(42) 
(fc',6) 
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and similarly for I{a,k'^). The impact factor H42|l comes 
from the region of the longitudinal momenta in the emit- 
ted gluon close to n. Note that the integral in the r.h.s 
of Eq. (|42|l is convergent and behaves like Xnk^lP' for 
k oo. Numerically, the correction to the impact 
factor H42|l is quite significant - it is responsible for the 
difference between the dotted and dash-dot lines in Fig. 
[TUl below. (The NLO corrections to the virtual photon 
impact factor calculated recently also appear to be big, 
see Refs. mil. 



C. Unpolarized dipole-dipole scattering 

The "unpolarized" dipole-dipole scattering corre- 
sponds to the cross section (p^ averaged over the ori- 
entation of the dipoles : 



a{a,b;s) = / ^' ^^^^ 



Note that in unpolarized 7*7* scattering the impact fac- 
tor lASp depends only on \a\ and |5| so the cross section 
of the unpolarized 7*7* scattering is given by the 
"unpolarized" dipole-dipole cross section in Eq. (|43|l in- 
tegrated over the dipole sizes with weights being impact 
factors 



cr(a, b; 7]) = /(iV^ - 1) coth^ t] 



(44) 




{k + k')^ 



27r 



d^kd^k' [1 - ./o(fca)][l - Mk'b)] 

l^[l~Mka)][l-Mkb)] 

f 2 k^ 
7] coth rj d p 



p'^{k — p)-^ 



1/5 2nf \ , 1 , 2 , 

— ^ In ^ H In n^ab 

47rV3 SNcJ F 27r ^ 



The explicit form of the function af(k, fc'; ry) is presented 
in the Appendix B. 



VI. NUMERICAL ESTIMATES OF THE 
DIPOLE-DIPOLE CROSS SECTION 



For numerical estimates, let us consider the unpolar- 
ized scattering of equal-size dipoles. In this case we have 
only one-scale problem, so it is natural to take ^ = 
where a is the dipole size. The amplitude H44|) reduces 
to 

(T(a, a; rj) = lena'^alia) coth^r/ 1-^60^(0) $(7/) (45) 



where 

= Stt / d'^kd'^k' [1 - Jo(fc)][i - Jo(fc')] 
'd^k 



(46) 



- Stt 



[i-Mk)Y 



rj coth r] 



/M af{k,k';T]) 
[k + k'f 

d^p Ink 



p'^{k-pY 27rfc2 



where we have rescaled k and k' by a and took nf — 
Nc — 3- This expression should be compared to the 
asymptotical form (|40|l which gives 



,{a,a;r]) = 167ras(o)o^ l + 6as{a)^ 



I'LLA+IF [ 



iv) (47) 



where 

„ fd^kd^k' , ^ ,,,, 
*LLA+iF = Stt / — — — [1 - Jo(fc)] 



X [1- Jo(fc')] 



fc2 k 

2t] \Ti{k + k' f /k"^ (2k'^ 



{k + k'f {k + k'f -k"^ \k 
ln(fc + k' f/k'^ (2k'^ - l^l - 2^^ ^'^ 



(48) 



-i- 



{k + k'f -k''-" V F 
'd^k 



k'^k'^ 



8n 



k^ 



Ink 



p'^{k — pf 2tt 



is a sum of the LLA term ~ 77 and the impact factor term 
^const. Numerically, $lla+if(?7) — 1.65(77 — 3.25) so 



'asy 



{a,a;ri) = IGna ag{a) 



1 



88as{a){T]-3.25) 

(49) 

where 77 corresponds to the LLA approximation and - 
3.251 comes from the "dipole impact factor" H42|l . We 
have not calculated the next ~ i term in the asymptotic 
expansion H49|l but the fit to the Figure [TOI suggests the 
coefficient of order of 100 in front of it. 

On Figure ^1 we have depicted the comparison be- 
tween the function ^{rj) describing the exact second- 
order cross section, its asymptotics $lla+if(77) Eq. 
and the "pure" BFKL asymptotics <i>LLA(?y) = 1-647/. 

We have used the VEGAS algorithm of the Monte- 
Carlo technique with 10% accuracy to calculate the func- 
tion ^{rj). 

We see that the BFKL asymptotics starts rather late, 
at 77 '--^ 5, which translates into y/s ~ 10 GeV for the 
scattering of dipoles with the p-meson size a ~ 0.25/777. 
This is perhaps not surprising since even within the 
(coUinearly improved (27,]) NLO BFKL approximation it- 
self the asymptotics starts relatively late, at 7/ 6 -f- 8 
for the realistic dipoles |2^. It should be emphasized, 
however, that these are two different theoretical ques- 
tions relevant to the BFKL description of experimental 
cross sections: (i) how good is the BFKL approximation 
to the true high-energy asymptotics and (ii) when this 
true asymptotics starts making sense. For the second- 
order dipole-dipole scattering the first question is when 
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FIG. 10: Exact cross section and its BFKL asymptotics. The 
solid line is Eq. 11461 , the dotted one is $lla+if Eq. (I48II and 
the dashed line is ^lla = 1.65r; 



the impact factor correction ~ const can be neglected 
(in comparison to the LLA Ins term) whereas the sec- 
ond question is when the ^ corrections disappear. In 
this paper, we mainly address the second question and 
our conclusion is that the true asymptotics (represented 
by $LLA+iF in this order) starts rather late. 



at ?7 ~ 5. Of course, it would be of more interest to 
check this statement for the scattering of virtual pho- 
tons rather than dipoles, but the corresponding diagram 
for the photon-photon scattering contains 4 loops so one 
should not expect this calculation to be performed in the 
nearest future. 

Our statement about the lateness of asymptotical be- 
havior refers to the scattering of equal-size dipoles. In 
many cases such as the deep inelastic scattering (DIS) 
one of the dipoles is small. The conventional assump- 
tion is that the propagation amplitude for the small-size 
dipole through a hadron is proportional to the "area of 
the dipole" multiplied by gluon parton density nor- 
malized at /i^ = 1/a^. This formula is valid for the 
asymptotically small dipoles, but it is unknown for how 
small dipoles and at what energies it starts making sense. 
This is an important question since in the analysis of 
the DIS from a nucleus based on the non-linear equa- 
tion for the dipole evolution 0, [s^] we assume that 
at the energy of a few GeV and size of order of satura- 
tion scale [aillllllllsl (2-^-3 GeV for LHC) the dipole 
matrix element between the nucleon states satisfies this 
approximate formula a^G{a^^). Since we cannot perform 
model-independent calculations of the dipole-hadron am- 
plitudes at present, we can consider the "target" being 
the second dipole of a size of a typical hadron. The am- 
plitude then is the a b limit of Eqn. (|^ . For this 
DIS from a color dipole we can answer the question of at 
what size of the spectator dipole the a^G approximation 
start making sense. The study is in progress. 
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VII. CONCLUSIONS AND OUTLOOK 

Our main result is that the scattering of dipoles is as 
suitable process to study the energy behavior of QCD 
amplitudes as, for instance, the scattering of virtual pho- 
tons or oniums. As we discussed above, the reason why 
this conclusion could have been different is the infinite 
length of Wilson lines forming the dipoles that leads to 
the mixture of the IR and UV divergencies. There are 
individual diagrams where the upper cutoff in the inte- 
grals over the transverse momenta is s rather than - 
. . . . " 

or the normalization scale /i. Such contributions would 

lead to the additional non-BFKL In s terms coming from 
the integrals over the transverse momenta. Fortunately 
in the sum of all diagrams, such terms cancel so that 
the transverse integrals are cut with the dipoles sizes or 
normalization point fj,. The resulting dipole-dipole am- 
plitude is an analytic function of the angle between the 
dipoles so that in principle, one can get the BFKL kernel 
from the Euclidean calculation. 

The second conclusion is that for the dipole-dipole 
scattering the high-energy asymptotics starts rather late. 
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APPENDIX A: VIRTUAL PHOTON IMPACT 
FACTOR 



In the momentum representation, the 7* impact factor 
for the transverse photons with the initial polarization ca 
and final is given by 



/-^(fcx) = ^Jdadp 



1 



{(l-2aa)(l-2^/3)(eA, 6^4)^1 



(Al) 



+ 4aa/3/3 



2{eA,k){e'^,k) - {eA,e'^)kl\'j 
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where x = 1 — x. If we average over the transverse po- 
larizations, the last term drops so that 



dad(3 



(1 - 2aa)(l - 2f3l3)kj 



(A2) 



In the coordinate representation, the impact factor (|A2|I 
takes the form 



/■^(.xx) ^ 
I^ix^) 

X 2Ki 



4^ 



dadf3 



(l-2aa)(l- 2/3/3) 
Anaa 



-Pa 



X 



,2 ^2 



9 2 



aa 

w 



„ , aa 



(A3) 



where X is the McDonald function. Note that the spin- 
averaged impact factor (|A3p depends only on x\. 



APPENDIX B: EXPLICIT FORM OF THE 
SECOND-ORDER AMPLITUDES 

As we mentioned above, the amplitude is a sum of five 
functions coming from the five terms in the product of 
two Lipatov vertices 



The explicit form of the first function is 



(Bl) 



(B2) 



2?7Cothry 



cosh r] 



AviTj cosh rj — {\ — ri — rj)^ 



2rj tanh rj \ Ari sinh 77 + 3 + — 



(1 — Ti) sinh 77 
2ry tanh rj \ Arj sinh^ r] + 3 + rj ~ ri 



( 1 — ) sinh 77 



Ar]{l + r, + rj) 
rirj sinh 2ri 

in coth ri ' 



1/2 



(l-r,0(l-r,) 



Iry^ - H- r,, -f rj 



where r, 



(fci+fe'i) 



In 



T and 

2„/„ I /„2 



-l-2sinh^,7(n + ^rf + ^ + 



* sinh'' jj 



(B3) 



The second function is 

A2 (B4) 
M(ri,r'i) - M(r2,r'i) - M(ri, r^) -f M(r2, r^) 



= -2rt 
where 



(ri - r2)[r\ ~ r'^) 



M(r,,r,) 



(B5) 



1 



^ (2ri cosh^ ?7 

4 cosh rjrirj — (1 — — rj)^ 

+ 1 — r.; — ^j^li -|- (2rj cosh^ 77 + 1 — — rj)/^ 



277 sinh 277 + iTT ( 2 cosh 77 



1 - ri - r.,- 



^ ) sinh?7 



and rt 



Note that it does not 



-t _ (fci~fc2)i 

contribute to the forward scattering and therefore to the 
total cross section. 

The third and the fourth functions are given by 



A. = 



4 tanh 77 



(ri - r2)(r'i - r^) 



(B6) 



(ri + r'i)Af(ri, ri) - (n + r^)M(ri,r^) 
(ra + r'^)M (ra, r^) - (ra + r^)Af (r2, r^) 



and 

^4 



(B7) 



2 tanh 77 
ri - r2 
2 tanh 77 



iV(r;,7'i) +iV(r^,ri) - (n ^ r2) 
N{r^,r\) + N{r2,r\)-{r\^r'2) 



where 

^{r,,rj] 



Arirj cosh 7/ — (1 — — rj)^ 



(B8) 



4ri cosh^ ?7 — rj 



{ri coth 77 — 277 + ztt) 



n -1 

-|- 277 — ITT + (1 + — rj)lj coth 7; 



/ Ticosh 77 o- 

dlTT-l / 5 ilTT. / C0Sh?7 

y Tj sinh 7/ + 1 Y 



Finally, 

A5 = - 



1/1 



rir2 vr';^ 
Zir2 - Z2?'i 



1 / 1 



r'jr2 vri r2 

'l'2 '2'l 



''i?'2(''i ^ ^2) sinh 7/ r'j^r2(r'j^ — Tj) sinh 7y 



(B9) 



All the functions Ai are analytical functions of the energy 
(= rj) . In the Euclidean region rj should be replaced with 
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iO where 9 is the angle between the dipoles (= Wilson 
rectangles). It is easy to see that the functions Ai become 
real after such a replacement. 

In the case of dipole-dipole cross section we need the 
real parts of the functions ljB2|l - (|B9|1 at fcj — and 
k\ = k' 2. After some algebra we get 

af{ry-ri) (BIO) 
= r'; 77) + /3(r, r'; ry) + Ji{r, r';f]) + f^{r, r';r]) 

Here /i(r, r'; 77) = ^A{r, r'; rf) and 

/3 = 2tanh2 7?^-^(r + r')m(r,r';?7) (Bll) 
— n(r ,r;r7) + — n(r,r ; 77) 



k 

h 



-2tanhr7 



dr' 



^ l-rdl/dr ^ ^ l-r'dl'/dr' r + r' 
2r^ sinh^ 77 2r'^ sinh^ r; 



where m{r,r';ri) 
^N{r, r';r]). 



^M{r^r';ri) and n{r,r']'q) 
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